\(\int \cot ^4(e+f x) \sqrt {a-a \sin ^2(e+f x)} \, dx\) [466]

Optimal result

Integrand size = 26, antiderivative size = 91 \[ \int \cot ^4(e+f x) \sqrt {a-a \sin ^2(e+f x)} \, dx=\frac {2 \sqrt {a \cos ^2(e+f x)} \csc (e+f x) \sec (e+f x)}{f}-\frac {\sqrt {a \cos ^2(e+f x)} \csc ^3(e+f x) \sec (e+f x)}{3 f}+\frac {\sqrt {a \cos ^2(e+f x)} \tan (e+f x)}{f} \]

[Out]

2*csc(f*x+e)*sec(f*x+e)*(a*cos(f*x+e)^2)^(1/2)/f-1/3*csc(f*x+e)^3*sec(f*x+e)*(a*cos(f*x+e)^2)^(1/2)/f+(a*cos(f
*x+e)^2)^(1/2)*tan(f*x+e)/f

Rubi [A] (verified)

Time = 0.15 (sec) , antiderivative size = 91, normalized size of antiderivative = 1.00, number of steps used = 5, number of rules used = 4, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.154, Rules used = {3255, 3286, 2670, 276} \[ \int \cot ^4(e+f x) \sqrt {a-a \sin ^2(e+f x)} \, dx=\frac {\tan (e+f x) \sqrt {a \cos ^2(e+f x)}}{f}-\frac {\csc ^3(e+f x) \sec (e+f x) \sqrt {a \cos ^2(e+f x)}}{3 f}+\frac {2 \csc (e+f x) \sec (e+f x) \sqrt {a \cos ^2(e+f x)}}{f} \]

[In]

Int[Cot[e + f*x]^4*Sqrt[a - a*Sin[e + f*x]^2],x]

[Out]

(2*Sqrt[a*Cos[e + f*x]^2]*Csc[e + f*x]*Sec[e + f*x])/f - (Sqrt[a*Cos[e + f*x]^2]*Csc[e + f*x]^3*Sec[e + f*x])/
(3*f) + (Sqrt[a*Cos[e + f*x]^2]*Tan[e + f*x])/f

Rule 276

Int[((c_.)*(x_))^(m_.)*((a_) + (b_.)*(x_)^(n_))^(p_.), x_Symbol] :> Int[ExpandIntegrand[(c*x)^m*(a + b*x^n)^p,
x], x] /; FreeQ[{a, b, c, m, n}, x] && IGtQ[p, 0]

Rule 2670

Int[sin[(e_.) + (f_.)*(x_)]^(m_.)*tan[(e_.) + (f_.)*(x_)]^(n_.), x_Symbol] :> Dist[-f^(-1), Subst[Int[(1 - x^2
)^((m + n - 1)/2)/x^n, x], x, Cos[e + f*x]], x] /; FreeQ[{e, f}, x] && IntegersQ[m, n, (m + n - 1)/2]

Rule 3255

Int[(u_.)*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]^2)^(p_), x_Symbol] :> Int[ActivateTrig[u*(a*cos[e + f*x]^2)^p]
, x] /; FreeQ[{a, b, e, f, p}, x] && EqQ[a + b, 0]

Rule 3286

Int[(u_.)*((b_.)*sin[(e_.) + (f_.)*(x_)]^(n_))^(p_), x_Symbol] :> With[{ff = FreeFactors[Sin[e + f*x], x]}, Di
st[(b*ff^n)^IntPart[p]*((b*Sin[e + f*x]^n)^FracPart[p]/(Sin[e + f*x]/ff)^(n*FracPart[p])), Int[ActivateTrig[u]
*(Sin[e + f*x]/ff)^(n*p), x], x]] /; FreeQ[{b, e, f, n, p}, x] &&  !IntegerQ[p] && IntegerQ[n] && (EqQ[u, 1] |
| MatchQ[u, ((d_.)*(trig_)[e + f*x])^(m_.) /; FreeQ[{d, m}, x] && MemberQ[{sin, cos, tan, cot, sec, csc}, trig
]])

Rubi steps \begin{align*} \text {integral}& = \int \sqrt {a \cos ^2(e+f x)} \cot ^4(e+f x) \, dx \\ & = \left (\sqrt {a \cos ^2(e+f x)} \sec (e+f x)\right ) \int \cos (e+f x) \cot ^4(e+f x) \, dx \\ & = -\frac {\left (\sqrt {a \cos ^2(e+f x)} \sec (e+f x)\right ) \text {Subst}\left (\int \frac {\left (1-x^2\right )^2}{x^4} \, dx,x,-\sin (e+f x)\right )}{f} \\ & = -\frac {\left (\sqrt {a \cos ^2(e+f x)} \sec (e+f x)\right ) \text {Subst}\left (\int \left (1+\frac {1}{x^4}-\frac {2}{x^2}\right ) \, dx,x,-\sin (e+f x)\right )}{f} \\ & = \frac {2 \sqrt {a \cos ^2(e+f x)} \csc (e+f x) \sec (e+f x)}{f}-\frac {\sqrt {a \cos ^2(e+f x)} \csc ^3(e+f x) \sec (e+f x)}{3 f}+\frac {\sqrt {a \cos ^2(e+f x)} \tan (e+f x)}{f} \\ \end{align*}

Mathematica [A] (verified)

Time = 0.08 (sec) , antiderivative size = 47, normalized size of antiderivative = 0.52 \[ \int \cot ^4(e+f x) \sqrt {a-a \sin ^2(e+f x)} \, dx=-\frac {\sqrt {a \cos ^2(e+f x)} \left (-3-6 \csc ^2(e+f x)+\csc ^4(e+f x)\right ) \tan (e+f x)}{3 f} \]

[In]

Integrate[Cot[e + f*x]^4*Sqrt[a - a*Sin[e + f*x]^2],x]

[Out]

-1/3*(Sqrt[a*Cos[e + f*x]^2]*(-3 - 6*Csc[e + f*x]^2 + Csc[e + f*x]^4)*Tan[e + f*x])/f

Maple [A] (verified)

Time = 0.79 (sec) , antiderivative size = 75, normalized size of antiderivative = 0.82

method	result	size
default	\(-\frac {\cos \left (f x +e \right ) a \left (3 \left (\sin ^{4}\left (f x +e \right )\right )+6 \left (\sin ^{2}\left (f x +e \right )\right )-1\right )}{3 \left (\cos \left (f x +e \right )-1\right ) \left (1+\cos \left (f x +e \right )\right ) \sin \left (f x +e \right ) \sqrt {a \left (\cos ^{2}\left (f x +e \right )\right )}\, f}\)	\(75\)
risch	\(\frac {i \sqrt {\left ({\mathrm e}^{2 i \left (f x +e \right )}+1\right )^{2} a \,{\mathrm e}^{-2 i \left (f x +e \right )}}\, \left (-3 \,{\mathrm e}^{8 i \left (f x +e \right )}+36 \,{\mathrm e}^{6 i \left (f x +e \right )}-50 \,{\mathrm e}^{4 i \left (f x +e \right )}+36 \,{\mathrm e}^{2 i \left (f x +e \right )}-3\right )}{6 \left ({\mathrm e}^{2 i \left (f x +e \right )}-1\right )^{3} f \left ({\mathrm e}^{2 i \left (f x +e \right )}+1\right )}\)	\(105\)

[In]

int(cot(f*x+e)^4*(a-a*sin(f*x+e)^2)^(1/2),x,method=_RETURNVERBOSE)

[Out]

-1/3*cos(f*x+e)*a*(3*sin(f*x+e)^4+6*sin(f*x+e)^2-1)/(cos(f*x+e)-1)/(1+cos(f*x+e))/sin(f*x+e)/(a*cos(f*x+e)^2)^
(1/2)/f

Fricas [A] (verification not implemented)

none

Time = 0.28 (sec) , antiderivative size = 66, normalized size of antiderivative = 0.73 \[ \int \cot ^4(e+f x) \sqrt {a-a \sin ^2(e+f x)} \, dx=-\frac {{\left (3 \, \cos \left (f x + e\right )^{4} - 12 \, \cos \left (f x + e\right )^{2} + 8\right )} \sqrt {a \cos \left (f x + e\right )^{2}}}{3 \, {\left (f \cos \left (f x + e\right )^{3} - f \cos \left (f x + e\right )\right )} \sin \left (f x + e\right )} \]

[In]

integrate(cot(f*x+e)^4*(a-a*sin(f*x+e)^2)^(1/2),x, algorithm="fricas")

[Out]

-1/3*(3*cos(f*x + e)^4 - 12*cos(f*x + e)^2 + 8)*sqrt(a*cos(f*x + e)^2)/((f*cos(f*x + e)^3 - f*cos(f*x + e))*si
n(f*x + e))

Sympy [F]

\[ \int \cot ^4(e+f x) \sqrt {a-a \sin ^2(e+f x)} \, dx=\int \sqrt {- a \left (\sin {\left (e + f x \right )} - 1\right ) \left (\sin {\left (e + f x \right )} + 1\right )} \cot ^{4}{\left (e + f x \right )}\, dx \]

[In]

integrate(cot(f*x+e)**4*(a-a*sin(f*x+e)**2)**(1/2),x)

[Out]

Integral(sqrt(-a*(sin(e + f*x) - 1)*(sin(e + f*x) + 1))*cot(e + f*x)**4, x)

Maxima [A] (verification not implemented)

none

Time = 0.33 (sec) , antiderivative size = 57, normalized size of antiderivative = 0.63 \[ \int \cot ^4(e+f x) \sqrt {a-a \sin ^2(e+f x)} \, dx=\frac {8 \, \sqrt {a} \tan \left (f x + e\right )^{4} + 4 \, \sqrt {a} \tan \left (f x + e\right )^{2} - \sqrt {a}}{3 \, \sqrt {\tan \left (f x + e\right )^{2} + 1} f \tan \left (f x + e\right )^{3}} \]

[In]

integrate(cot(f*x+e)^4*(a-a*sin(f*x+e)^2)^(1/2),x, algorithm="maxima")

[Out]

1/3*(8*sqrt(a)*tan(f*x + e)^4 + 4*sqrt(a)*tan(f*x + e)^2 - sqrt(a))/(sqrt(tan(f*x + e)^2 + 1)*f*tan(f*x + e)^3
)

Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 3443 vs. \(2 (83) = 166\).

Time = 1.82 (sec) , antiderivative size = 3443, normalized size of antiderivative = 37.84 \[ \int \cot ^4(e+f x) \sqrt {a-a \sin ^2(e+f x)} \, dx=\text {Too large to display} \]

[In]

integrate(cot(f*x+e)^4*(a-a*sin(f*x+e)^2)^(1/2),x, algorithm="giac")

[Out]

-1/24*(48*(sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/
2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x)
*tan(1/2*e)^2 - sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*t
an(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2
*f*x) - 2*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2
*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*e))/(
(tan(1/2*f*x)^2 + 1)*(tan(1/2*e)^2 + 1)) + (3*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1
/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/
2*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x)^5*tan(1/2*e)^10 + 3*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*
x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4
- 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x)^4*tan(1/2*e)^11 + sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*ta
n(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1
/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x)^3*tan(1/2*e)^12 - 18*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*
e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e
)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x)^5*tan(1/2*e)^8 - 51*sqrt(a)*sgn(tan(1/2*f*x)^
4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)
*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x)^4*tan(1/2*e)^9 - 30*sqrt(a)*sgn(tan
(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*ta
n(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x)^3*tan(1/2*e)^10 - 3*sqrt(
a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2
*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x)^2*tan(1/2*e)^11
+ 72*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x
)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x)^4*ta
n(1/2*e)^7 + 177*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*
tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/
2*f*x)^3*tan(1/2*e)^8 + 93*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f
*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) +
1)*tan(1/2*f*x)^2*tan(1/2*e)^9 + 3*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 -
tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan
(1/2*e) + 1)*tan(1/2*f*x)*tan(1/2*e)^10 + 18*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/
2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2
*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x)^5*tan(1/2*e)^4 + 72*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x
)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 -
4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x)^4*tan(1/2*e)^5 - 186*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*
tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan
(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x)^2*tan(1/2*e)^7 - 114*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/
2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2
*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x)*tan(1/2*e)^8 - 2*sqrt(a)*sgn(tan(1/2*f*x)^4
*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*
tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*e)^9 - 3*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1
/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/
2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x)^5*tan(1/2*e)^2 - 51*sqrt(a)*sgn(tan(1/2*f*
x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f
*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x)^4*tan(1/2*e)^3 - 177*sqrt(a)*sgn
(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) -
4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x)^3*tan(1/2*e)^4 - 186*
sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*ta
n(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x)^2*tan(1/2*
e)^5 + 42*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2
*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*e)^7
+ 3*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^
3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x)^4*tan(
1/2*e) + 30*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4 - 4*tan(1
/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*tan(1/2*f*x
)^3*tan(1/2*e)^2 + 93*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2*f*x)^4
- 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e) + 1)*t
an(1/2*f*x)^2*tan(1/2*e)^3 + 114*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan
(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/
2*e) + 1)*tan(1/2*f*x)*tan(1/2*e)^4 + 42*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)
^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x
)*tan(1/2*e) + 1)*tan(1/2*e)^5 - sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan
(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/
2*e) + 1)*tan(1/2*f*x)^3 - 3*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 - tan(1/2
*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan(1/2*e)
+ 1)*tan(1/2*f*x)^2*tan(1/2*e) - 3*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e)^3 -
tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f*x)*tan
(1/2*e) + 1)*tan(1/2*f*x)*tan(1/2*e)^2 - 2*sqrt(a)*sgn(tan(1/2*f*x)^4*tan(1/2*e)^4 - 4*tan(1/2*f*x)^3*tan(1/2*
e)^3 - tan(1/2*f*x)^4 - 4*tan(1/2*f*x)^3*tan(1/2*e) - 4*tan(1/2*f*x)*tan(1/2*e)^3 - tan(1/2*e)^4 - 4*tan(1/2*f
*x)*tan(1/2*e) + 1)*tan(1/2*e)^3)/((tan(1/2*f*x)^2*tan(1/2*e) + tan(1/2*f*x)*tan(1/2*e)^2 - tan(1/2*f*x) - tan
(1/2*e))^3*tan(1/2*e)^3))/f

Mupad [B] (verification not implemented)

Time = 17.07 (sec) , antiderivative size = 364, normalized size of antiderivative = 4.00 \[ \int \cot ^4(e+f x) \sqrt {a-a \sin ^2(e+f x)} \, dx=\frac {\left (\frac {1{}\mathrm {i}}{f}-\frac {{\mathrm {e}}^{e\,2{}\mathrm {i}+f\,x\,2{}\mathrm {i}}\,1{}\mathrm {i}}{f}\right )\,\sqrt {a-a\,{\left (\frac {{\mathrm {e}}^{-e\,1{}\mathrm {i}-f\,x\,1{}\mathrm {i}}\,1{}\mathrm {i}}{2}-\frac {{\mathrm {e}}^{e\,1{}\mathrm {i}+f\,x\,1{}\mathrm {i}}\,1{}\mathrm {i}}{2}\right )}^2}}{{\mathrm {e}}^{e\,2{}\mathrm {i}+f\,x\,2{}\mathrm {i}}+1}+\frac {{\mathrm {e}}^{e\,3{}\mathrm {i}+f\,x\,3{}\mathrm {i}}\,\sqrt {a-a\,{\left (\frac {{\mathrm {e}}^{-e\,1{}\mathrm {i}-f\,x\,1{}\mathrm {i}}\,1{}\mathrm {i}}{2}-\frac {{\mathrm {e}}^{e\,1{}\mathrm {i}+f\,x\,1{}\mathrm {i}}\,1{}\mathrm {i}}{2}\right )}^2}\,8{}\mathrm {i}}{f\,\left ({\mathrm {e}}^{e\,2{}\mathrm {i}+f\,x\,2{}\mathrm {i}}-1\right )\,\left ({\mathrm {e}}^{e\,1{}\mathrm {i}+f\,x\,1{}\mathrm {i}}+{\mathrm {e}}^{e\,3{}\mathrm {i}+f\,x\,3{}\mathrm {i}}\right )}+\frac {{\mathrm {e}}^{e\,3{}\mathrm {i}+f\,x\,3{}\mathrm {i}}\,\sqrt {a-a\,{\left (\frac {{\mathrm {e}}^{-e\,1{}\mathrm {i}-f\,x\,1{}\mathrm {i}}\,1{}\mathrm {i}}{2}-\frac {{\mathrm {e}}^{e\,1{}\mathrm {i}+f\,x\,1{}\mathrm {i}}\,1{}\mathrm {i}}{2}\right )}^2}\,16{}\mathrm {i}}{3\,f\,{\left ({\mathrm {e}}^{e\,2{}\mathrm {i}+f\,x\,2{}\mathrm {i}}-1\right )}^2\,\left ({\mathrm {e}}^{e\,1{}\mathrm {i}+f\,x\,1{}\mathrm {i}}+{\mathrm {e}}^{e\,3{}\mathrm {i}+f\,x\,3{}\mathrm {i}}\right )}+\frac {{\mathrm {e}}^{e\,3{}\mathrm {i}+f\,x\,3{}\mathrm {i}}\,\sqrt {a-a\,{\left (\frac {{\mathrm {e}}^{-e\,1{}\mathrm {i}-f\,x\,1{}\mathrm {i}}\,1{}\mathrm {i}}{2}-\frac {{\mathrm {e}}^{e\,1{}\mathrm {i}+f\,x\,1{}\mathrm {i}}\,1{}\mathrm {i}}{2}\right )}^2}\,16{}\mathrm {i}}{3\,f\,{\left ({\mathrm {e}}^{e\,2{}\mathrm {i}+f\,x\,2{}\mathrm {i}}-1\right )}^3\,\left ({\mathrm {e}}^{e\,1{}\mathrm {i}+f\,x\,1{}\mathrm {i}}+{\mathrm {e}}^{e\,3{}\mathrm {i}+f\,x\,3{}\mathrm {i}}\right )} \]

[In]

int(cot(e + f*x)^4*(a - a*sin(e + f*x)^2)^(1/2),x)

[Out]

((1i/f - (exp(e*2i + f*x*2i)*1i)/f)*(a - a*((exp(- e*1i - f*x*1i)*1i)/2 - (exp(e*1i + f*x*1i)*1i)/2)^2)^(1/2))
/(exp(e*2i + f*x*2i) + 1) + (exp(e*3i + f*x*3i)*(a - a*((exp(- e*1i - f*x*1i)*1i)/2 - (exp(e*1i + f*x*1i)*1i)/
2)^2)^(1/2)*8i)/(f*(exp(e*2i + f*x*2i) - 1)*(exp(e*1i + f*x*1i) + exp(e*3i + f*x*3i))) + (exp(e*3i + f*x*3i)*(
a - a*((exp(- e*1i - f*x*1i)*1i)/2 - (exp(e*1i + f*x*1i)*1i)/2)^2)^(1/2)*16i)/(3*f*(exp(e*2i + f*x*2i) - 1)^2*
(exp(e*1i + f*x*1i) + exp(e*3i + f*x*3i))) + (exp(e*3i + f*x*3i)*(a - a*((exp(- e*1i - f*x*1i)*1i)/2 - (exp(e*
1i + f*x*1i)*1i)/2)^2)^(1/2)*16i)/(3*f*(exp(e*2i + f*x*2i) - 1)^3*(exp(e*1i + f*x*1i) + exp(e*3i + f*x*3i)))